Self-energy effects in electronic Raman spectra of doped cuprates due to magnetic 

fluctuations 



m 
o 

(N 



^1" 



C3 
O 
O 

Oh- 
33 ■ 



C3 
O 

o 



> 
in 

ON 

in 

rn 
O 
m 



x 



Roland Zeyher 1 and Andres Greco 2 

1 Max- Planck- Institut fur Festkorperforschung, Heisenbergstrasse 1, D-70569 Stuttgart, Germany 
2 Departamento de Fisica, Facultad de Ciencias Exactas e Ingenieria and IFIR(UNR-CONICET), 2000 Rosario, Argentina 

(Dated: March 26, 2013) 

We present results for magnetic excitations in doped copper oxides using the random phase 
approximation and itinerant electrons. In the [1,0] direction the observed excitations resemble 
dispersive quasi-particles both in the normal and superconducting state similar as in recent resonant 
inelastic X-ray scattering (RIXS) experiments. In the [1,1] direction the excitations form, except 
for the critical region near the antiferromagnetic wave vector Q = (ji, n), only very broad continua. 
Using the obtained spin propagators we calculate electron self-energies and their effects on electronic 
Raman spectra. We show that the recently observed additional peak at about twice the pair breaking 
in Bi ff symmetry below T c in HgBa2Cu04+5 can be explained as a self-energy effect where a broken 
Cooper pair and a magnetic excitation appear as final states. The absence of this peak in B2 9 
symmetry, which probes mainly electrons near the nodal direction, is explained by their small self- 
energies compared to those in the antinodal direction. 

PACS numbers: 74.25. nd,74.25. Ha, 74.72.Gh 



Electronic Raman spectra in doped cuprates are dom- 
inated by a pair breaking peak in the superconducting 
state whereas in the normal state they are rather struc- 
tureless. The main features of the spectra in the su- 
perconducting state can be explained within BCS theory 
using free quasi-particles. [l[ On the other hand RIXS 
experiments have shown the existence of spin fluctua- 
tions in doped cuprates in form of dispersive, broad but 
well-defined excitations, both in the normal and the su- 
perconducting state. 0, Q The interaction between spin 
fluctuations and electrons gives rise to self-energies which 
have been used to explain, for instance, the observed 
kinks in the electronic dispersion in many high-T c su- 
perconductors. 0-@| Self-energy effects should also be 
present in electronic Raman spectra causing for instance 
the recently observed broad peak above the pair break- 
ing peak in HgBa 2 Cu04 +( 5. 7] It is the purpose of this 
Letter to calculate the properties of spin fluctuations, the 
resulting self-energies and their effects on Raman scatter- 
ingand to compare the results with recent experiments. 
@J3| Our calculations are based on an itinerant picture 
of magnetism and will be limited to the optimally doped 
and overdoped region to avoid unsolved problems with 
the existence and nature of the pseudogap in the under- 
doped regime. 

The Raman susceptibility can be written in the super- 
conducting state as, Q 

Xa(w m ) = ^ 7a(k)(IIll,u(k, iv m ) + II 22 , 22 (k, jv m ) 



-IIi2,2i(k,ii/ m ) - n 2 i,i2(k,if m )), (1) 



/OC 
dede'Aij(k, e)A k i(k, 
-OC 



iv m + e' — e ' 
(2) 



The indices assume the values 1 or 2, v m de- 

notes a bosonic Matsubara frequency, / the Fermi func- 
tion, and N c is the number of primitive cells. As- 
suming a square lattice and using the lattice constant 
as the length unit 7 Q (k) stands for the form factors 
7i(k) = t(cos k x + cos k y ) + At' cos k x cos k y , 73 (k) = 
t(cos(k x — cos/cy), 74 (k) = — At' sink x sinfc y , correspond- 
ing to the representations Ai g , Bi 9 , and B 2g of the point 
group T>4h , respectively, t and i! are effective nearest and 
second-nearest neighbor hopping amplitudes. 

Using the Nambu representation the inverse of the 2x2 
electron Green's function matrix G[k,iw n ) is 

G _1 (k,iw n ) = (iu) n - S+(k, iw n ))r 
-(e(k)+E_(k,tw n ))7S-A(k)n. (3) 
Ti and T3 are Pauli matrices, To the 2x2 unit matrix, and 
oj n a fermionic Matsubara frequency. £±(k, W n ) stands 
for (Eu(k,ia>„) T E u (-k, —iw n ))/2, where E u (k,io; n ) 
is the (1,1) component of the 2x2 self-energy matrix. e(k) 
is the electron energy 

e(k) = — 2i(cos k x + cos k y ) — At' cos k x cos k y — /1, (4) 



and the spectral function Ay is defined by 
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iff). The self-energy E is calcu- 



lated from the Fock diagram with a bosonic propagator 
describing spin fluctuations. To simplify the calculation 
we do not consider the non-diagonal part of E assuming 
that it has already been taken into account by the 
phenomenological gap parameter A(k). Assuming also 
that the frequency-independent Fock term has been 
accounted for in the phenomenological hoppings t and t' 
the frequency-dependent part of En is given by [4| 

Eu(k, iu n ) = 

3 (0), 



q, iw„ 



0- (5) 
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x(k, iv m ) is the zz component of the spin susceptibility 
which assumes in the RPA the form 0, 



X(k, iv„ 



X (0) (k^ m ) 



1- J(k) X (°)(k,ii/ m )- 



(6) 



X^°^(q, w m ) represents the bare spin bubble. J(k) is 
equal to 2J(cos 

kx ~t~ cos fcy ) where is the Heisenberg 
constant. Since self-energy effects in Raman spectra are 
expected to be rather small we may use instead of 
Gn in Eq.©. 

In the following calculations we used the parameters 
t' /t = —0.25 and /i/t = —0.9 corresponding to the dop- 
ing S = 0.15, where t is the effective nearest-neighbor 
hopping amplitude and about 0.15eV in optimally doped 
YBCO. We expect a similar value in the case of 
HgBa2Cu04+5- J was determined so that the experimen- 
tal value of 0.056 eV 11| for the magnetic resonance en- 
ergy in the superconducting state was reproduced yield- 
ing J/t — 1.30. The bare order parameter A was fixed 
to A/t = 0.45 which reproduces roughly the position of 
the observed pair breaking peak in Bi ff symmetry in the 
sample with the highest T c of 94K of Ref. Q. 

An important ingredient for the calculation of £ is the 
dynamical spin susceptibility x(k, w). Its imaginary part 
"(k,uj) has been discussed extensively in the literature 
in conjunction with the resonant mode, i.e., in the su- 
perconducting state and near the antiferromagnetic wave 
vector Q. We extend these calculations to wave vectors 
throughout the Brillouin zone (BZ) and also to the nor- 
mal state. Such calculations are interesting in view of 
recent RIXS measurements which cover a more ex- 
tended region. The upper and lower diagrams in Fig. 1 
show x"(k, along the [1, 1] direction in the normal and 
the superconducting state, respectively. For wave vectors 
away from Q (and also from the origin (0,0)) the spectra 
are extremely broad and structureless, exhibiting often 
several flat maxima. Only the center of gravity of the 
lines indicates an increase with momentum up to around 
the middle of the zone and then a pronounced decrease 
towards Q. In the normal state the curve for momentum 
Q shows a well-pronounced critical peak due to the prox- 
imity of the antiferromagnetic instability. In the super- 
conducting state this peak moves to somewhat higher en- 
ergies and becomes the resonant peak with a finite width 
due to our employed damping constant r\ = 0.1. Being a 
bound state inside the gap it attracts a huge amount of 
spectral weight. This holds to a lesser degree even for the 
momentum 3Q/4 where the structureless continuum in 
the normal state transforms into a well-pronounced peak 
in the presence of superconductivity. 

x"(k, lu) looks quite different in the [1,0] direction as 
illustrated in Fig. 2. It shows both in the normal and 
the superconducting state well-defined and narrow peaks 
which disperse as a function of momentum similar as the 
spin waves of a Heisenberg model. However, in our case 
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FIG. 1. (Color online) x"(k,i^>) along the [1,1] direction in 
the normal (upper diagram) and superconducting (lower dia- 
gram) state at T=0. 



the sharp peaks do not describe collective excitations. 
This means that the sharp peaks in Fig. 2 do not corre- 
spond to poles of the denominator of x(k, lu) but are due 
to dynamical nesting properties of x^ (k, w) . The energy 
scale of the magnetic dispersion is in our case determined 
by the electron band dispersion, i.e., by the kinetic en- 
ergy of the electrons and not by the Heisenberg constant 
J. This is especially clear at the point (it, 0): In the RPA 
J drops out in x and the energy of the magnetic exitation 
is soley due to the kinetic energy of the electrons. Using 
different band parametrizations (for instance, those given 
in Table I of Ref. [l(| ) we found very similar results as in 
Figs. 1-2 showing that these results are generic for doped 
cuprates. 

x"(k, ui) has recently been determined by RIXS mea- 
surements along the [1, 0] direction for several copper ox- 
ides and successfully interpreted within a local-moment 
model including damping. Q Fig. 2 and its comparison 
with Ref. Q shows that the nature of the magnetic ex- 
citations, their dispersion and their band width are sim- 
ilar in the itinerant and in the local moment approach, 
though the underlying phyics is quite different. One rea- 
son for this is that J and the effective hopping t have 
very similar values of about 150 meV. In the [1, 1] di- 
rection the two approaches may show larger differences 
because the itinerant model predicts in this case rather 
broad and irregular line shapes except for a small critical 
region around (7r,7r). 

Fig. 3 depicts the imaginary part of £(k, cj) (upper 
diagram) and the spectral function ^4n(k, uj) (lower di- 
agram) for the momentum on the Fermi line along the 
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FIG. 2. (Color online) along the [1,0] direction in 

the normal (upper diagram) and superconducting (lower dia- 
gram) state at T=0. 
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FIG. 3. Imaginary part of the self-energy (upper diagram) 
and spectral function (lower diagram) for the momentum on 
the Fermi line in [1, 1] direction in the normal and supercon- 
ducting state. 



[1,1] direction. In the normal state Jm,S(k,w) is zero 
at w = 0, increases away from zero due to scattering 
of electrons with spin fluctuations, and finally decreases 
due to the finite band width. In the superconducting 
state electron scattering is impeded by the gap leading 
to a depiction of spectral weight around u) = in im£ . 
On the other hand superconductivity causes a large shift 
of spectral weight of spin fluctuations towards low ener- 
gies, especially near (tt, tt). Momentum transfers of about 
(tt, 7r) lead in general to scattered electrons far away from 
the Fermi line for the considered intial electron momen- 
tum. The result is the superconductivity-induced peak 
around u/t ~ 2., i.e., well away from the Fermi energy. 
The spectral function Au(k,ui) consists of a simple peak 
at uj = without any sidebands. 

Fig. 4 shows the same as Fig. 3 but for the electron 
momentum on the Fermi line along the (tt, 0) — (tt, tt) di- 
rection. Now there are important electronic transitions 
between the surroundings of (tt, 0) and (0, tt) with mo- 
mentum transfers near (tt, tt), i.e., where spin fluctuations 
have a large spectral weight and low energies. In the nor- 
mal state the peak due to critial scattering near (n, tt) in 
the upper part of Fig. 1 causes the low-lying peaks in 
the imaginary part of E(k,Lj), as seen in Fig. 4. In the 
superconducting state the same scattering processes are 
even much stronger due to the appearance of the resonant 
mode which carries a large amount of spectral weight. As 
a result the distance between the two low-lying peaks in 
im£(k,Li;) in the normal state increases due to the su- 



perconducting gap and at the same time the peaks gain 
a lot of strength. The corresponding spectral functions 
reflect the low- lying peaks in form of sidebands, which 
are strongly pronounced especially in the superconduct- 
ing state. The position of the sideband peaks would be 
in the absence of dispersion the sum of half of the su- 
perconducting gap plus the resonance frequency. The 
dispersion of electrons and of the resonance peak causes 
an additional shift to the positions shown in the figure. 

Finally, Fig. 5 gives our results for the Raman inten- 
sity —Imxa(w + iv) where a denotes one of the three 
symmetries Bi g , B 2s and Ai g . In each case the dashed 
line corresponds to the normal and the solid line to the 
superconducting state. Let us first discuss the Bi g spec- 
trum (upper part in Fig. 5). In the free case, defined 
by £ — 0, the Raman intensity would consist in the nor- 
mal state of a <5-function at u> = (not shown in the 
figure) and a Lorentzian-like pair breaking peak at the 
unrenormalized energy 2A. Turning on E(k, w) the <5- 
function is replaced by a continuum in the normal state 
which increases rather abruptly from zero at zero fre- 
quency to a finite value and then decreases very slowly 
towards higher frequencies. This is the result of elec- 
tron scattering from spin fluctuations. The solid curve, 
referring to the superconducting state and Bi 9 symme- 
try, exhibits two maxima. The first one is due to the 
renormalized pair breaking peak, the second one due to 
the additional emission of a magnetic excitation, for in- 
stance, the resonance mode. The position of the addi- 
tional high-frequency peak can also be estimated from 
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FIG. 4. Imaginary part of the self-energy (upper diagram) 
and spectral function (lower diagram) for the momentum on 
the Fermi line in (tt, 0) — (it, ir) direction in the normal and 
superconducting state. 
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FIG. 5. (Color online) Raman spectra with Bi 9 (upper part), 
B29 and Ai s symmetries (lower part) in the normal (dashed 
lines) and superconducting (solid lines) state. 



the spectral function An. According to Eq.([T]) the Ra- 
man intensity at frequency ui is essentially be obtained 
as a folding of two spectal functions displaced by the fre- 
quency w. As a result the pair breaking peak arises from 
tansitions between the two dominating peaks in Fig. 4. A 
smaller peak is expected from transitions between one of 
the large peaks and a sideband peak yielding a frequency 
of about 1.3 in good agreement with Fig. 5. The ratio of 
the positions of the two peaks is about 1.9 which is near 



the experimental value. 0] We stress that the position 
and the intensity of the second peak in this spectrum is 
not fitted but a consequence of the parameters J and A 
determined such that the observed resonance frequency 
and the pair breaking peak are at least approximately re- 
produced. Using the electronic parameters of Ref. [4{ the 
additional peak in the B lff spectrum is practically absent. 
The reason for this is that the self-energy in this case is 
substantially smaller than in our case. The difference in 
the magnitude of the self-energies can at least partially 
be traced back to the large difference in the resonance fre- 
quencies, namely 56 meV in HgBa 2 Cu04 + 5 and 40 meV 
in YBCO for nearly optimally doped samples. 

The B-2 9 spectrum (contained in the lower diagram in 
Fig. 5) is in the normal state similar to that of the B lg 
symmetry except that its intensity is smaller by one or- 
der of magnitude. In the superconducting state it shows 
a pair breaking peak but no additional peak at higher 
frequencies, in agreement with experiment. 12[ The rea- 
son for this can easily be understood: The form factor in 
the B2 9 symmetry weights electrons near the nodal di- 
rection strongly. There the gap and also the self-energy 
are small. As a result the pair breaking moves somewhat 
to lower frequencies compared to the Bi s symmetry and 
no additional peak appears as a result of self-energy ef- 
fects. For completeness we also show the A lg spectrum 
(lower diagram in Fig. 5) where Coulomb screening has 
been taken into account. Here we encounter the well- 
known and hitherto unsolved problem that its intensity 
is too small by about one order of magnitude compared 
to experiment. [13;] Thus we will exclude it from our 
discussion. 

We used in our treatment a broken Cooper pair times 
a spin fluctuation as final states. Using the same in- 
teractions one may also consider a competing process 
where two spin fluctuations or paramagnons appear as 
final states. For zero or very small doping such a 
localized description is certainly appropriate whereas in 
the optimal or overdoped regime a description in terms 
of itinerant electrons should be a better choice. For in- 
stance, it allows easily to understand why the additional 
high-energy peak shows a strong dependence on T c sim- 
ilar like the resonance mode [7] or why its existence is 
tied to the existence of the pair breaking peak and the 
depletion region below it. It also has been argued [l4| 
that two-magnon processes are neglegible small in the 
Bi s and B2 ff configurations and thus may be discarded 
in these cases. 

In summary, we have shown that an itinerant picture 
of magnetism may account for the observed dispersive 
magnetic excitations in [1, 0] direction in doped cuprates 
and also for the observed additional structure above the 
pair breaking peak in electronic Raman scattering in 
HgBa 2 Cu0 4+ 5. 
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